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Abstract 

Dense rubbery networks are highly entangled polymer systems, 
with significant topological restrictions for the mobility of neighbour- 
ing chains and crosslinks preventing the reptation constraint release. 
In a mean field approach, entanglements are treated within the fa- 
mous reptation approach, since they effectively confine each individual 
chain in a tube-like geometry. We apply the classical ideas of repta- 
tion dynamics to calculate the effective rubber-elastic free energy of 
anisotropic networks, nematic liquid crystal elastomers, and present 
the first theory of entanglements for such a material. 
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1 Introduction 



Rubbery polymer networks are complex randomly disordered amor- 
phous systems. The simplest theoretical models consider them as 
being made of "phantom chains", where each chain is thought to be 
a three-dimensional random walk in space. To form a network, the 
chains are crosslinked to each other at their end points, but do not 
interact otherwise, in particular they are able to fluctuate freely be- 
tween crosslinks. This has an unphysical consequence that the strands 
can pass through each other. If one tries to avoid this assumption, the 
theory is confronted with the intractable complexity of entanglements 
and their topological constraints. The mean field treatment of entan- 
gled polymer melts and semi-dilute solutions is the classical reptation 
theory Jl], ||] going back to the early seventies, which has been a spec- 
tacular success in describing a large variety of different physical effects. 
However, the parallel description of crosslinked rubbery networks has 
been much less successful. First of all, there is a significant differ- 
ence in entanglement topology: in a melt the confining chain has to 
be long enough to form a topological knot around a chosen polymer; 
even then the constraint is only dynamical and can be released by a 
reptation diffusion along the chain path. In a crosslinked network, 
any loop around a chosen strand becomes an entanglement, which 
could be mobile but cannot be released altogether. A number of other 
complexities arise from the coupling between imposed deformations 
and chain anisotropy, the stress-optical effects [||, f|, || and nematic 
interactions between chain segments |], [jj. 

In addition, the polymer network can be spontaneously anisotropic, 
forming a liquid crystalline elastomer (LCE). This area has attracted 
a significant experimental and theoretical interest in recent years. In 
nematic LCE, the strands preferably orient themselves along one di- 
rection, forming a nematic liquid crystal order. The response to an 
external deformation is now of a much richer nature, with antisymmet- 
ric stress and internal torques depending on the relative angle of the 
director to the axis of deformation (8|. Liquid crystalline elastomers 
combine remarkable properties of both its components, liquid crys- 
tals and rubbers, but also show physical properties that place them 
in a separate category from any other material. Several new physical 
phenomena have been discovered in LCE: (a) spontaneous, reversible 
shape changes of up to 400% on temperature change; (b) "soft elas- 
ticity" - mechanical deformation, involving modifications of internal 
nematic microstructure, without (or with very low) stress; (c) me- 
chanical and electric instabilities involving director reorientation, in 
special cases discontinuous jumps; (d) solid phase nematohydrody- 
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namics and unusual rheology, leading to anomalous dissipation and 



acoustic effects. Recent reviews [|9|, [10|] describe the current state of 
affairs in this field. Our challenge in this paper is to bring the micro- 
scopic theoretical description of nematic rubbers on the same level as 
in the classical isotropic rubbers, in particular, to account for chain 
entanglements. 

An early model of elastic response of entangled rubbers was de- 
veloped by Edwards [11]: in tradition with the melt theory, it as- 
sumed that the presence of neighbouring strands in a dense network 
effectively confines a particular polymer strand to a tube, whose axis 
defines the primitive path. Within this tube, the polymer is free to ex- 
plore all possible configurations, performing random excursions, par- 
allel and perpendicular to the axis of the tube. One can show that on 
deformation of the sample the length of the primitive path increases. 
Since the arc length of the polymer is constant, the amount of chain 
available for perpendicular excursions is reduced, leading to a reduc- 
tion in entropy and an increase in rubber-elastic free energy. A number 
of further attempts have been made to derive a self-consistent theory 
of entangled rubber elasticity. Of this list, the most significant are the 
scaling "localisation model" of Gaylord and Douglas ll]| , the "slip-link 
model" of Ball, Doi, Edwards and Warner (BDEW) p3|, accounting 
for entanglements as local mobile confinement sites linking two inter- 
wound strands, and the "hoop model" by Higgs and Ball (HB) |I~4} ], 
who assumed that entanglements localise certain chain segments. In 
an article developing the reptation theory of rubber elasticity for clas- 
sical isotropic networks |l5| , we give a more detailed overview of these 
and other theories. 

In our current work, we extend the tube model to treat the elas- 
ticity of anisotropic networks of liquid crystalline polymers. To our 
knowledge, this is the first time that a reptation model has been 
applied to treat the effects both of the entanglements and of the 
anisotropic nature of the nematic network. The tube model provides 
a more accurate microscopic description in the sense that it keeps 
track of the allocation of chain segments and their excursions between 
the points of entanglement. The next Section briefly reviews the ideal 
phantom-network approach to the elasticity of nematic rubber and in- 
troduces the tube model and its properties, giving the full expression 
for nematic rubber-elastic free energy. Section 3 contains the discus- 
sion of the model and its results, including the linear-response limit. 
We conclude by comparing the results of the present theory with those 
of the ideal phantom network and analyse which physical properties 
of LCE seem to be most sensitive to the effect of chain entanglements. 
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2 Nematic elastomer network 



Before developing our model for macroscopic elasticity of densely en- 
tangled rubber, we briefly review the well-known results of the phan- 
tom chain network theory, which provides the basics to most other 
theoretical models. 

Phantom chain approximation 

Assuming that a single polymer performs a free random walk in three 
dimensions, one finds that the end-to-end distance Rq obeys a Gaus- 
sian distribution in the long chain limit. This result goes back far in 
history: one can review its derivation and consequences in the classical 
text on this subject H. The distribution of Ro is given by 

where b is the monomer step length and N the number of steps of one 
chain trajectory. 

In a nematic polymer, irrespective of the particular mesogenic 
mechanism, the monomer steps acquire a preferred orientation along 
the director no- Accordingly, the end-to-end distance distribution 
function of a strand becomes anisotropic as well: 

/ q \ 3/2 

p <^> - (saw) < det ^~ 1/2 w 

x<*p(-^»o- ST' ■*<>)■ 

where Nb = L is the contour length of the chain, and the matrix 
Uptakes account of the anisotropy: 

(| )ij = lofaj + {Iq ~ lo)n 0i n j. 

This matrix of anisotropic chain steps is directly measurable from the 
average chain shape, given by (RiRj) = \{§g)ijL. The principal values 

of this effective step-length matrix, Iq and ijj, reflect the spontaneous 
nematic order in the material. In the isotropic phase, e.g., above the 
nematic transition temperature T n \, ll = Iq = b and one trivially 
recovers the isotropic Gaussian distribution ([[]). The difference Iq — 
Iq is proportional to the nematic order parameter Q. However, the 
explicit form of this dependence is different in different models of 
nematic polymers. In a most simple case of freely jointed chain of 
rods of length b one obtains Iq = 6(1 + 2Q), Iq = b(l — Q), while 
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in the hairpin regime of semiflexible main-chain nematic polymer the 
anisotropy could become very large: /jj oc exp[3/(l — Q)], cf. Jl6| ]. The 
power of the ideal theory of nematic rubber elasticity S is in that 
it is independent of such model considerations and only uses a single 
model parameter - the ratio r = Iq/Iq, or equivalently, r = / {R±) 
for the principal values of the gyration radius. We shall see below that 
this attractive feature is reproduced in the theory of entangled nematic 
networks. 

The entropic free energy of such an anisotropic random walk is 
given by the logarithm of the number of conformations with the fixed 
Rq and has the form 

(3F = — InP(Ro) = — flj ■ | Q 1 • Ro + const, 

where f3 = 1/ksT is the inverse Boltzmann temperature. At forma- 
tion of the network, i.e. at crosslinking, the polymer melt is assumed 
to obey the anisotropic Gaussian distribution (|2|), which is then per- 
manently frozen in the network topology. In the phantom network ap- 
proximation, the lateral restrictions on the chain thermal motion are 
neglected and different network strands interact only at the crosslink- 
ing points. 

One then assumes that the junction points deform affinely with 
respect to their initial positions Rq by the macroscopic deformation 
A, hence we can write R = XRq. Therefore, the deformation A alters 
the free energy of each strand. The change of free energy per chain of 
the whole network can be calculated by the usual quenched averaging: 

PF = -(lnP(Ai2o))p ( flo) = ^Tr(| Q • A T -I" 1 • A), (3) 

where we have dropped an irrelevant constant and found a new ex- 
pression for the chain step-length matrix after deformation: 

with the rotated director, (n • no) = cos 6, and possibly changed prin- 
cipal values and l ± . The overall elastic free energy density, in the 
first approximation, is simply (]3|) multiplied by the number of elasti- 
cally active network strands in the system n c h per unit volume, which 
is proportional to the crosslinking density: 

F el = ^Tr(| .A T -|^-A), (4) 
with the rubber modulus fi = n^ksT, cf. || for detail. 



5 



The phantom-network model of rubber elasticity is a popular first 
approximation. There are several reasons for its overall success in 
spite of obvious oversimplifications. The crosslinking points connect 
the ends of different strands together and thus reduce local fluctuations 
- and, therefore, alter the single chain statistics. However, in spite of 
an apparent complexity of this problem, it has been shown [17] that 
this effect merely introduces a trivial multiplicative factor of the form 
1 — 2/(j), where <f> is the junction point functionality. Secondly, one 
can assume that the deformation preserves the volume, since the bulk 
(compression) modulus is by a factor of at least 10 4 greater than the 
shear modulus, which is proportional to fi; this implies the constraint 
det A = 1. Thirdly, the quenched average in equation (|3|) does not 
average over chains of different arc lengths, but the fact that the result 
is independent of arc length, generalises the result to apply for chains 
of arbitrary length, or even for a polydisperse ensemble of chains. In 
the particular case of nematic LCE, this simple model provides a rich 
crop of theoretical predictions described in greater detail in quoted 
review articles. 



The tube model 

Following the original ideas of Edwards [jllj] , we assume that one par- 
ticular network strand is limited in its lateral fluctuations by the pres- 
ence of neighbouring chains. Therefore each segment of a given strand 
only explores configurations in a limited volume, which is much smaller 
than in the random coil state. Hence, the whole strand fluctuates 
around a mean path, which we call the primitive path. Effectively, 
the chain is confined to exercise its thermal motion only within a tube 
around the primitive path due to the presence of neighbouring chains. 
This primitive path itself can be considered as a random walk with an 
associated typical step length, which is much bigger than the monomer 
step length |jT8[ . The step length of the primitive path divides the tube 
into tube segments, as sketched in Fig. |], and therefore determines the 
number of tube segments M along one polymer strand. 

Note that all chains are in constant motion, altering the local con- 
straints they impose on each other. Hence, the tube is a gross sim- 
plification of real situation. However, one expects this to be an even 
better approximation in rubber than in a corresponding melt (where 
the success of reptation theory is undeniable), because the restriction 
on chain reptation diffusion in a crosslinked network eliminates the 
possibility of constraint release. 

To handle the tube constraint mathematically, we assume that the 
chain segments are subjected to a quadratic potential, restricting their 
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Figure 1: A nematic polymer strand is surrounded by neighbouring chains, 
which effectively confine the strand to a tube. The tube segment m contains 
s m monomer steps, where the index m runs from 1 to M. Since the direction 
of chain steps is, on average, uniaxial (as illustrated by the drawing of an 
ellipsoidal shape of a free chain gyration on the right), the number of steps 
in each tube segment depends on the orientation of A m with respect to the 
nematic director n. 

motion transversely to the primitive path. Along one polymer strand 
consisting of N monomers of effective step length b, there are M tube 
segments, each containing s m , m = 1,...,M monomer steps. We 
infer the obvious condition 

M 

Y J S m = N. (5) 

m=l 

In effect, one has two random walks: the topologically fixed primitive 
path and the polymer chain restricted to move around it - both having 
the same end-to-end vector Rq, between the connected crosslinking 
points. 

Each tube segment m can be described by the span vector A m , 
joining the equilibrium positions of the strand monomers at the two 
ends of each tube segments. The number of tube segments M (or, 
equivalently, the number of chain entanglements, M — 1) is a free 
parameter of the theory, ultimately determined by the length of each 
polymer strand and the average "entanglement density" . 

Since the primitive path is a topologically frozen characteristic of 
each network strand, we shall assume that all primitive path spans 
A m deform affmely with the macroscopic strain: A^ = AA m . This 
is the central point in the model: the rubber elastic response will 
arise due to the change in the number of polymer configurations in a 
distorted primitive path. To evaluate the number of conformations, 
we look separately at chain excursions parallel and perpendicular to 
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the tube axis, within each span A m . Effectively, this amounts to 
introducing a new coordinate system for each tube segment, with one 
preferred axis along A m . The constraints exerted by the other chains 
only constrain the considered polymer in its lateral motion. Hence, 
we recover the behaviour of a one-dimensional random walk in the 
direction of A m , giving rise to Gaussian statistics in the long chain 
limit. Note that only one third of the steps s m in this tube is involved 
in the longitudinal excursions. We therefore obtain for the number of 
longitudinal excursions in a tube segment m, cf. Fig. [l], 

The spontaneous anisotropy of nematic polymer chain is reflected in 
@ by accounting for the difference in number of chain conformations 
in a given tube segment, depending on its orientation with respect to 
the local nematic director (the principal axis of step- length matrix 

To determine the number of transverse excursions, we introduce 
the Green's function for the steps made by the chain in the plane 
perpendicular to the local tube axis A m . In effect, we consider a two- 
dimensional random walk, with a total number of steps (2s m /3), in a 
centrosymmetric quadratic potential. For each of these two perpen- 
dicular coordinates, the Green's function satisfies the following modi- 
fied diffusion equation (see e.g. and its extension for the uniaxial 
nematic case in The argument that follows, about transverse 

chain movements in a confining potential, has a very simple conclu- 
sion - that the conformational effects are irrelevant for calculation of 
rubber elasticity and only important factor is the number of chain 
steps attributed to this degree of freedom. Equally, the anisotropic 
(nematic) nature of chain random walk does not contribute to the 
entropy of strongly confined transverse excursion. Although the full 
anisotropic treatment is possible, here we shall use a much shorter 
and transparent version of the isotropic chain confined in the tube; its 
Green's function satisfies the differential equation for each of the two 
coordinates: 

d 1 d 2 1 \ 

ds ~ 2 &2 ft? + 2 q ° x2f J G ( Xi > x f> s ) = S ( x f ~ Xi ) 5 ( s )> ( 7 ) 

where the Xi and x f are the initial and final coordinates of the random 
walk with respect to the tube axis and qo determines the strength 
of the confining potential. The equation (|^) is very common in the 
physics of polymers and its exact solution is known. However, we only 
need to consider a particular limit qobs m ^> 1 of this solution, which 
is the case of dense entanglements (resulting in a strong confining 
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potential) and/or of a large number s m of monomers confined in the 
tube segment. Outside this limit, that is, when the tube diameter 
is the same order as the arc length of the confined chain, the whole 
concept of chain entanglements becomes irrelevant. In the strongly 
confined limit the solution has a particularly simple form [11 1: 



G(xi, xf, s) oc exp ( -^{xf + xj) - ^q bs m ) . (8) 



Remembering that there are two coordinates describing the transverse 
excursions, we obtain for the two-dimensional Green's function of the 
tube segment m: 

G m {n, Tf, s) oc exp (-^qobsm) exp (-^( r i + > ( 9 ) 

where r$ and rj are the initial and final transverse two-dimensional 
coordinates. 

The total number of transverse excursions is proportional to the 
integrated Green's function: 



W^ T) oc / dri J dr f G m {r u r f ;s) 



Since the Green's function (^) does not couple the initial or final co- 
ordinates to the number of segments s m , this integration will only 
produce a constant normalisation factor which can be discarded. Ex- 
actly the same conclusion is reached if the modified diffusion equation 
for anisotropic chain is considered. 

Gathering the expressions for statistical weights of parallel and 
perpendicular excursions (and returning to the fully anisotropic de- 
scription) , one obtains the total number of configurations of a polymer 
segment consisting of s m monomers in a tube segment of span A m : 

W m = Wg>wW (10) 

Therefore, we find for the full number of configurations of the whole 
strand 

W = j dsf- [ ds M (l[ W m ) S ( s m - n\ (11) 

q q \m=l / \m=l / 

where we have implemented the polymer contour length constraint 
The statistical summation in (|Tl|) takes into account the reptation 
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motion of the polymer between its two crosslinked ends, by which the 
number of segments, s m , constrained within each tube segment can be 
changed and, thus, equilibrates for a given conformation of primitive 
path. 

Rewriting the delta- function as 5(x) = ^ / dk e lkx , we proceed by 
finding the saddle points which make the exponent of statistical 
sum stationary. It can be verified that the normalisation factors 
l/y/sm contribute only as a small correction to the saddle points 



k 2ba qQ b + \k) 



(12) 



The integral in ( |Tl| ) is consequently approximated by the steepest 
descent method. We repeat the same procedure for the integration 
of the single auxiliary variable k, responsible for the conservation of 
the polymer arc length. The saddle point value k* , inserted into fll2|), 
gives 

A^o 1/2 A m | 

Sm ~ T M \r 1 ' 2 A I ' 

which is the equilibrium number of steps the nematic polymer makes 
in a tube segment characterised by the axis vector A m . By complet- 
ing the saddle-point integration, we finally obtain the total number 
of configurations of one strand, confined within a tube whose prim- 
itive path is described by the set of vectors {A m }. The statistical 
weight W associated with this state is proportional to the probability 
distribution: 

W(A 1 ,...,A M )<xP({A m }) (13) 

2 



■X 



exp (-2w(Em=il^ 1/2 A m |) -hobN 



Af-1 



The scalar \t^ 1 / 2 A m \ reflects the length of the m-th step of the prim- 
itive path, modified by its projection on the uniaxial matrix of chain 
step- lengths. This expression is a result parallel to the ideal Gaussian 
P(Rq) in equation (||) for a unentangled chain. Note that the chain 
end-to-end distance Rq is also the end-to-end distance of the primitive 
path random walk: J2m=i A m = Ro- 

Free energy of deformations 



From the equation (|13|) we obtain the formal expression for free energy 
of a chain confined to a tube with the primitive path conformation 
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{A m }, 0F = -hxW, or 




+(M 



where we have dropped irrelevant constants arising from normalisa- 
tion. We now perform a procedure which is analogous to the one used 
to obtain equation (^) . In the polymer melt before crosslinking, we as- 
sume that the ensemble of chains obeys the distribution in (13) giving 
the free energy per strand fllj). The process of crosslinking not only 
quenches the end points of each of the crosslinked strands, but also 
quenches the nodes of the primitive path A m , since the crosslinked 
chains cannot disentangle due to the fixed topology of the network. In 
our mean field approach, the tube segments described by A m are con- 
served. For evaluating the quenched average, note that the statistical 
weight ( |i~3| ) treats all tube segments m in an equivalent way. This 
allows one to perform the summation over the index m, separating 
the diagonal and the off-diagonal terms: 



+M(M -l)(|| 1 /2 Am |||- 1 /2 Ar 
+(M -1) An f£lio~ 1/2A -l)y 



for arbitrary values of m and n ^ m; the brackets (• • •) refer to the 
average with the probability P({A m }) given in (13). 

Any mechanical deformation expressed by the general strain tensor 
A will affinely transform A m into A m = AA m . It could also affect 
the nematic order: the director n could adopt a different orientation 
under deformation and the degree of average chain anisotropy r may 
change as well. In other words, the matrix which characterises 
the anisotropy of the steps, transforms into a new matrix i„ with 
different eigenvalues and I 1 - in a reference frame rotated by the 
angle 9. Hence |£~ 1//2 A m | transforms into \i e X ^ 2 AA m | on deformation, 
but the distribution P({A m }) remains unchanged. Bearing this in 
mind, we can evaluate the averages (|l5|), leading to the free energy 
per crosslinked chain. The Appendix gives a more detailed account of 
how one evaluates the averages. The resulting elastic energy density 
takes the form 

Fel = ^ ^J-1 T A T i 

3 P 3M + 1 v=0 = = 9 - v ; 



11 



+ M (m-i)ih|^ i/2 -a-i; /2 i, 

where we use the notations: 

\S 1/2 -h-^\ = h /^ 1/2 -A-|f e| (17) 



|e]=l 

Inli^-A-if I = ^ / dfi In II- 1 / 2 • A • if e| (18) 

|e|=l 

(the overline notation 7TT refers to the angular averaging over the orien- 
tations of an arbitrary unit vector e used to contract a corresponding 
matrix into a vector, before calculating its absolute value). 

Expressions (17) and ([TJ) can be evaluated in various particular 
cases of deformation A an d director orientation. The Appendix gives 
a result for uniaxial deformation along the director, where A takes a 
diagonal form with A" = A and X ± = l/y/\. Explicit formulae for ( p"7D 
and (|l^) need to be inserted into (|fj) to give the full elastic energy. 



3 Discussion 



From the expression (16), we can recover the elastic free energy of an 
ideal phantom-chain nematic network by taking the case M = 1. This 
limit means physically that the polymer strand is placed in one single 
tube, tightly confined to the axis. Mathematically, a random walk in 
three dimensions with iV steps is equivalent to a random walk in one 
dimension along a given direction with N/3 steps. This fact is the 
underlying reason why we recover the phantom chain network result 
by taking M = 1 in our model. 

On the other hand, as the number of tube segments M becomes 
large, one obtains a rubber-elastic elastic energy of the form 

F A = ixM (((||, 1/2 • A • i 1/2 |» 2 + (In || e - 1/2 • A • i 1/2 |>) • (19) 

There are two ways to have a physical situation corresponding to this 
limit of M > 1: either the polymer melt is very dense, causing a high 
entanglement density, or the polymer chain is very long between its 
crosslinked ends. In the latter case, the polymer strand experiences 
many confining entanglements along its path. 

Recall that the F e \ is the elastic energy density, which relates to 
the free energy per chain F as: F e \ = n c h_F, where n c h is the density 



12 



of crosslinked strands. We can assume that in a polymer melt, the 
chain density is inversely proportional to the volume of an average 
chain, hence inversely proportional to the contour length of this chain: 
n c h oc 1/L. In case of the phantom chain network the rubber modulus 
H = n^ksT [equation (0)]. One concludes in this case that the elastic 
energy F e \ scales with 1/L, and therefore -F e i — > as the chains become 
infinitely long! This unphysical behaviour reflects the fact that the 
phantom chain model assumes the entanglement interactions of the 
chains irrelevant. Clearly, this assumptions breaks down in the long 
chain limit, where one expects the entanglements to play a crucial 
role. 

This unphysical behaviour is overcome by our expression (|l^). As 
the strands become longer, they will experience more entanglements, 
generating more confining tube segments. We could reasonably as- 
sume that the number of entanglements and therefore the number of 
tube segments scales linearly with the strand length L: M oc L. Con- 
sidering expression (|l9|) , one can note that the corresponding rubber 
modulus does not vanish in the limit L — » do, but remains a constant 
corresponding to the "rubber plateau" in a densely entangled melt. 

Considering the particular case of uniaxial strain along the con- 
stant nematic director n, one can examine one of the key physical ef- 
fects found in nematic elastomers - the spontaneous mechanical defor- 
mations as the degree of anisotropy is changed, for instance, by chang- 
ing the temperature (and thus the nematic order parameter Q(T) and 
the effective chain anisotropy r). Within the ideal phantom-chain 
model (||), applying a uniaxial deformation along the director n with 
A" = A and A^ = 1/y/X, one obtains 

where zjj and Zq~ are the principal values of £ , and similarly Z" and 
I 1 - the ones of I , the anisotropy of a state after the deformation (of 
course, in this case no director rotation occurs). The free energy is 
minimised by the strain 

A m = (/"/ ± //j|/ ± ) 1 / 3 , 

which describes a spontaneous uniaxial deformation of a nematic rub- 
ber, first discovered theoretically in [ff^j and mentioned in the lit- 
erature ever since. For instance, if the initial state ^ Q is isotropic (at 
T > T n i), then A m = (Z" /Z -1 -) 1 / 3 , a function of nematic order parameter 
Q(T) and could reach a remarkable value of 400% uniaxial extension 
in a highly anisotropic main-chain nematic rubber ||19|]. 



13 



This result is not altered by the complicated additional terms in 
( |j~6| ) : remarkably, exactly the same deformation A m minimises all three 
corresponding expressions derived from ([if]), which are given in the 
Appendix. 

If we now assume that both the chain anisotropy and the director 
n are kept fixed under the deformation, £g = i^, and that the strain 
tensor A is diagonal in the reference frame of the anisotropy matrix, 
then we observe that the matrices in (p]|)-(|i~8|) are all diagonal. Hence 
the anisotropic terms cancel each other out, and we are left with the 
same elastic energy as in the isotropic case Hence, even if the 

material is anisotropic, its linear elastic modulus does not depend on 
the orientation under the above assumptions of unchanged degree of 
anisotropy r = l^/l^: the Young's moduli E" = E 1 - = E. However, 
the modulus for simple shear is partially affected by the anisotropy of 
the nematic rubber. Consider Am = 5ij + euiVj, with u and v, the two 
orthogonal unit vectors defining the simple shear. If the deformation 
does not mix the parallel and perpendicular directions, i.e., if u and 
v are both perpendicular to n, then the shear modulus is the same as 
in the isotropic case, 

1 /4 2M + 1 I, , 11M + 5N , . 

G = -E = a ( — + -(M - 1 — - . (20) 

3 P \3 3M + 1 5 V ' 3M + 1 J y ' 

On the other hand, if one of the vectors u or v is parallel to the 
director n, then the shear modulus is changed by a factor of (Iq /Iq) 1 ^ 2 
or (Zq / Zq- ) 1 / 2 , respectively. 

If the material is not allowed to deform, any rotation of the nematic 
director n away from its equilibrium orientation no will cost energy. In 
phantom-chain networks, the trace formula (Q) gives the corresponding 
elastic free energy increase as a function of 9, the angle between n and 
n : 

(in the limit of small director rotation 6). This gives the expression 
for the relative rotation coefficient Di, first written down phenomeno- 



logically by de Gennes [2C] and extensively discussed in the literature 
[], [10[] . In the small strain limit A = 5 + e, the coupling between 
the director rotation u = [n x 5n] and the antisymmetric part of the 
strain J7j = eijkSjk can be written as: 

^Di [nx(fl- u)} 2 + D 2 n- g {s) ■ [n x (O - w)] , 

where is the symmetric part of the small strain. The entangle- 
ment model does not change the dependence D(r) qualitatively, but 
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introduces a coefficient associated with the entanglement density: 



Another key physical property of nematic rubbers is the effect of 
soft elasticity. Fundamental internal symmetries of an elastic medium 
with an independently mobile orientational degree of freedom, the 
nematic director n, demand that there is a particular relationship 
between the two relative rotation coefficients D\ and D2 and one of 
the linear shear moduli, C5, It has been shown [22] that there 



is a continuous set of such soft deformations (not necessarily small 
in amplitude), which by appropriately combining strains and director 
rotations can make the elastic response vanish completely: 

4 oft = lf^-lo- 1/2 . 

where U_ is an arbitrary unitary (3D rotation) matrix. It is quite 
obvious that substituting this strain tensor into the modified tube- 
model expression ( |l6|) will leave this free energy at its ground state 
level as well. It is, in fact, gratifying that these two crucial physical 
effects (thermal expansion and soft elasticity), which have attracted 
so much theoretical and experimental attention in recent years, are 
left intact within a much more complex theoretical description of a 
highly entangled nematic elastomer. 



Conclusion 

In this present work, we have analysed the behaviour of a uniaxial ne- 
matic polymer network in the presence of chain entanglements, which 
are treated within a tube model approach. We found that this leads 
to a significantly modified rubber-elastic energy which, in principle, 
should supersede the earlier molecular theory (||). The present model 
captures the physics of entanglements in a consistent way and, for the 
first time, takes into account an orientational effect of chain conforma- 
tion in the tube segments aligned at an arbitrary angle with respect 
to the uniform nematic director n. Since the role of entanglements 
is, from all points of view, much more significant in a crosslinked net- 
work, the theory provides a firmer ground for description of many 
theoretically known and experimentally tested results. 

We have to remark that our model only describes the equilibrium 
response of a network to deformation. Shortly after applying the defor- 
mation, the network will need to find a new microscopic equilibrium. 
Each polymer strand would redistribute the monomers between the 
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affinely modified tube segments, attributing more monomers to some 
segments, less to others, and eventually reaching a new optimal con- 
formation {s m }. This gives the expression for the rubber elastic free 
energy density (|i~6|). The dynamics of this relaxation is based on the 
sliding (reptation) motion along the primitive path while constraining 
the end points of it. This process would be reflected in a time de- 
pendence of the variable s m , which is the number of monomer steps 
attributed to the tube segment m. By describing this relaxation pro- 
cess, one could extend the present equilibrium model to describe the 
stress relaxation and the short time viscoelastic response of a nematic 
rubber. 

We appreciate many useful discussions with S.F. Edwards and M. 
Warner. S.K. gratefully acknowledges support from an Overseas Re- 
search Scholarship, from the Cambridge Overseas Trust and from Cor- 
pus Christi College. 

A Evaluation of quenched averages 

To evaluate averages (A^ • | X • A m ), (||~ 1/2 A m | ||~ 1/2 A n |) and 
(ln(^ |^q A m |)) in equation (|l5|), for the arbitrary m,n = 1,...,N, 
one needs to integrate the corresponding scalar functions of A m 
with respect to the probability distribution For this purpose, 

one has first to find the normalisation J\f of the distribution, which 
can most easily be achieved by introducing a new scalar variable 
u = Em=i l^o 1//2 Am| to simplify the exponent. It is also useful to 
change the integration variables from A m to a transformed vector 
A m = £~ 1//2 A m . One obtains then: 




u M-|Ai| 

J dAiA 2 J dA 2 A 2 • • • 
o o 
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-...-|Ai 



f-2| 

dA 



(«-E 



m= 



In the last step, we introduced spherical coordinates for the vari- 
ables A m , implemented the delta- function constraint u = E and 
used the fact that the variables A m are bound to be positive. The 
underlined expression is a multiple integral over the hyper-triangular 
domain in the space of {A m } and is a function of u, which we call 
Im{u). Since the integrals only involve power functions, Im(u) itself 
is a power in u. It is then evaluated via the iterative procedure, which 
generates the recursive relation and returns an explicit function: 

M-l 



,3Af-l 



f. 



with f M = n 



3m(3m + l)(3m + 2) 



The first two terms in the Eq. (|T5|) involve the diagonal (~ A^J or the 
off-diagonal (~ A m A n ) factors. In both cases, the integration proce- 
dure is analogous to that of normalisation factor M above, except that 
either one (m) or two (m 7^ n) integrals in the sequence contain an 
extra scalar factor of A m . The corresponding angular integration over 
the orientations of A m (producing a factor of 4ir in Jv) now becomes 
non-trivial, depending on its angle relative to tensors A and £ , when 
the sample is deformed. This angular integration is left unfinished 
here, since it depends on particular deformation and director geome- 
try; the main thermodynamic average of the diagonal (square) term 
returns the ideal trace- formula in the final free energy density (|l6|), 
while the off-diagonal average returns the expression (17). 
For the logarithmic term in (|l(^), one obtains: 




AA r 



Detg'/ 2 
N 



dA 1 A( 



du 



2b 2 N 



dA 2 A 2 2 ■ 



HA, 



f-ai 

dAw-i A 



u 



E 



M- 
m= 



n / <M r 

\m=l J 



In 



A X II" 1 / 2 - A ■t /2 e 1 \ + 



••• + A M _ 1 ||- 1/2 .A-^ /2 e 



Af-l 



E 



M- 

■m = 



)|^ 1/2 -A-^| 
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Here dVt m is the angular measure of orientations of the corresponding 
unit vector e m , along the modified tube segment vector A m . In the 
last term under the logarithm, the absolute value of A.m is substituted 
by its value from the delta- function constraint. The next step is to 
approximate the logarithm with its complicated angular-dependent 
argument: 

ln[A 1 ||- 1/2 -A-iy 2 Si| + ... 



+a m _ 1 |!, i/2 -a-i; /2 %/-i| 

+(n-^|A m |)|^ 1 / 2 -A-4 /2 eM|] 



InM + In ||~ 1/2 • A • |o /2 e M | + In 



1 + 



M-l 



+ E 



m=l 



1=0 ^ ' = '=(/ ®m\ 

||- 1/2 .A-|fe M | 



A 



rn 
U 



small since «>A m 
~ In \lg 1/2 ■ A • ll /2 e M \ + const. 

After this, all of the results of multiple integrals over dA m and dVt m 
cancel with the normalisation factor and the only relevant contribution 
arises from the angular integration of the scalar logarithmic term over 
the orientations of e^f, cf. equation (18). 

In the particular case of uniaxial deformation along the nematic 
director no, with A" = A (along no=const) and A^ = l/yf\, the 
evaluation of equations (|l6|)-(|T8|) gives: 



Tr(iL • A T 



A) 



\£e 1/2 -h-ll /2 \ 




In II, 



-1/2 



•a-4 /2 i 



ln(A) - 1 



arctan(^/y / r) 




with the notations 



Vx 3 



and 



If we assume that the anisotropy is not changed by the deformation, 
i.e. £ Q = £g, (and the director preserves its original orientation no) 
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then the elastic response of the nematic rubber is not different from 
isotropic behaviour [15]. 
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